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1 Introduction 

A synthesis of strictly twodimensional crystals composed of carbon atoms 
(graphene) [1] is promising a wealth of new phenomena and possible applica- 
tions in technology and industry [2]. The observation of anomalous transport 
properties, and, most exciting, the recent discovery of substantial field effect and 
magnetism at room temperature allows one to envisage graphene as a reasonable 
replacement of nanotubes in electronic applications, see, e.g., Refs. PIUG)]. 

By symmetry, the valence and conduction bands in graphene touch at the 
corners of the hexagonal Brillouin zone. In the vicinity of these points, the dis- 
persion relation is isotropic and linear, and the density of states at the Fermi level 
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is strictly zero, rising linearly in energy. An effective long-wavelength description 
of these electronic states can be written in terms of a continuum model which 
is based on the Dirac- Weyl equation for massless electrons in 2+1-dimensional 
space-time [6], [7] . 

Lowdimensional quantum systems of Dirac fermions can possess rather un- 
usual properties and, since the discovery of the effect of charge fractionalization 
[8], are generating current interest. Planar Dirac fermions in the background 
of the Abrikosov-Nielsen-Olesen vortex [HI [10] were studied in Ref. [11], and, 
recently, the results of this work have been used to consider the influence of 
the Kekule distortion in the honeycomb lattice on the electronic properties of 
graphene [121 13] • The present paper deals with yet another aspect, and our 
purpose is to elucidate the role of topological defects in the graphene lattice. 

Topological defects appear as a result of removing (inserting) one or several 
carbon atoms from (into) the honeycomb lattice without affecting the threefold 
coordination of other atoms; appropriately, the lattice surface is warped owing 
to positive (negative) curvature induced at the location of a defect. Assuming 
that the size of a defect is small as compared to the whole size of the graphene 
sample, our interest will be in the study of the influence of such a defect on the 
electronic properties of graphene. The consideration is based on the continuum 
model for long-wavelength electronic excitations, and various types of defects are 
characterized by just the number of carbon atoms removed or inserted; actually, 
the size of a defect is neglected. The graphene sheet with a defect takes shape 
of a cone with the value of the apex angle related to the number of removed 
atoms. When a defect is encircled, then sublattices, as well as inequivalent Fermi 
points, are entwined or left untwined, depending on the type of a defect. This 
imposes a condition on the electronic wave function on the graphene sheet: a 
phase is acquired under a rotation around a defect. If the phase commutes with 
the hamiltonian, then it can be eliminated by a singular gauge transformation 
which, on the other hand, introduces a fictitious point vortex that may be denoted 
in the following as a pseudomagnetic one. The flux of the pseudomagnetic vortex 
is related to the deficit (proficit) angle of the conical surface, i.e. to the number 
of removed (inserted) atoms. For certain types of defects the vortex flux takes 
fractional values in the units of 2tt. 

A general theory of planar relativistic fermionic systems in the background 
of a point magnetic vortex with arbitrary flux was proposed in Refs. [T^ [T5] ; 
in particular, the case of massless fermions was considered in Refs. [TBI [T7] . If 
the vortex flux is fractional in the units of the London flux, then an essentially 
irregular mode appears among the eigenmodes of the one-particle hamiltonian, 
and this in its turn gives rize to the appearance of an additional parameter - 
the self-adjoint extension parameter which specifies the boundary condition at 
the vortex point [TS|. The theory allows us to predict the density of states (not 
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local but total) [19] and the ground state quantum numbers, as well as their 
local densities [T5l [HI E] • The aim of the present paper is to apply this theory 
to graphene with a topological defect. The previous attempts to consider the 
electronic properties of graphene with topological defects in the framework of the 
continuum model approach [201 EH E2] have led to contradictory and, therefore, 
unconvincing results. In our opinion, this is due to the two circumstances: an 
inadequate choice of the condition for the electronic wave function in the case 
of the entwinement of sublattices and an inappropriate treatment of irregular 
eigenmodes of the one-particle hamiltonian. These deficiencies will be remedied 
in the present study. 

In the next section we review the derivation of the continuum model for pla- 
nar graphene in order to specify the notations used. In Section 3 we introduce 
topological defects in graphene in the framework of the continuum model, derive 
the condition for the electronic wave function in the case of an arbitrary defect, 
and demonstrate that the defect can be presented as a pseudomagnetic vortex 
at the apex of a graphitic nanocone. In Section 4 we consider the solution of 
the Dirac-Weyl equation for electronic excitations on a graphene sheet with a 
topological defect and find the density of states and the ground state charge. 
The results are discussed in Section 5. In Appendix A we show that the case of 
the three-pentagon defect coincides actually with that of the absence of defects. 
In Appendix B the method of a self-adjoint extension is applied to derive the 
condition for the irregular solution to the Dirac-Weyl equation. 

2 Continuum model for long- wave length 
electronic excitations 

Carbon atoms in graphene form a honeycomb lattice. The Bravais lattice is 
triangular, and the primitive cell is rhombic and contains two carbon atoms. If 
one atom is placed at the origin of the primitive cell, another one is displaced 
at d = (— d, 0), where d is the lattice spacing. Thus, the honeycomb lattice is 
composed of two triangular sublattices: sublattice Aa (black points) is generated 
by vectors = rijCi + m;C 2 , and sublattice A# (blank points) is generated by 

vectors = n^Ci + mjC 2 + d, where Ci = ^dj and c 2 = (^d, —^d^j are 

the basis vectors of the primitive cell (see Fig. la), and n i; m ; e Z (Z is the set 
of integer numbers). 

Each carbon atom in graphene has four valence electrons. As a result of 
hybridization, three of them build a-orbitals along the lattice links providing 
for their rigidity, while the fourth one makes 7r-orbital which is orthogonal to 
the lattice plane and is responsible for the conductive properties of graphene. 
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Retaining only nearest neighbour interactions in the tight-binding approximation, 
quantum-mechanical hopping of an electron on a planar honeycomb lattice is 
described with the use of hamiltonian 

3 3 

H = -tJ2 XV (rO&fc + uj) - t J2 XVfaMri + vj), (1) 

«eA A j=l ieA B j=l 

where t is the hopping amplitude, a^(rj) and a(r;) (fct(rj) and 6(rj)) are the 
creation and destruction operators acting on sublattice A^(A B ), which obey an- 
ticommutation relations 

[a(Ti), a t ( r i')]+ = [ & ( r i)> &t ( r i')]+ = $u', 

Uj(vj) are the triad of vectors which are directed to the nearest neighbours of an 
atom belonging to sublattice Aa(A b ), see Fig. lb, 



ui = (-d,0), u 2 



U3 — ( ^rf, — ^-d 



vi = (d,0), v 2 = , V 3=(-id,^d 



Using the Fourier transforms of the sublattice operators, 



a(Ti) = 



d 2 k 



e lkr *a(k), b(Ti) 



d 2 k 



ikr; 



&(k), 



(2) 



where Y is the first Brillouin zone, hamiltonian (1) is presented as 

d 2 k 



n 
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^(k)tf^(k), 



(3) 



where 



and 



^(k) = (a(k), 6(k)) T , V f (k) = (at(k), &t( k ) 
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Solving the eigenvalue problem, Hip(k) = Eip(k), one gets 
E 



\ j=i j'=i 



e ikv i' 



±t 



1 + 4 cos 



cos ( -k x d ) + cos I ^^k,,d 



(4) 



(5) 



As follows from Eq.(5), the one-particle energy spectrum consists of two surfaces 
(E > and E < 0) which intersect (E = 0) at six conical points that are located 
at 

k x = 0, k y = ±4n(3V3d)-\ 

k x = ±2n(3d)- 1 , k y = ±2tt(3 v / 3)^ 1 . 1 ' 

With one electron per site, the negative-energy states (valence band) are filled 
and the positive-energy states (conduction band) are empty, so the band structure 
is at half-filling with the Fermi level E = corresponding to six isolated points 
given in Eq.(6). The first Brillouin zone is a hexagon with corners identified with 
these Fermi points; among six of them, only two ones which can be taken as 
oppositely located are inequivalent, see Fig.2. 

Actually, H (4) has the meaning of the one-particle hamiltonian in the mo- 
mentum representation. The low-energy excitations can be studied by taking 
the continuum limit (d — > 0) at any of two inequivalent Fermi points. Choosing 
the pair of inequivalent points as K± = (0, ±47r(3v / 3d)~ 1 ) and keeping terms of 
order k — K±, one gets : 

H± = limrf- 1 #| k=K±+ft = !* ( _iK x °± Ky ^ o ^ ) = M-* 2 ** ± ^j/); 

(7) 

where v = ^th 1 is the Fermi velocity, and a 1 and a 2 are the off-diagonal Pauli 
matrices. Combining the contributions from K + and KL, one gets 

= hv [a X K x + a 2 Ky) , (8) 




where 

Making the Fourier transformation of Eq.(8), one gets the long- wavelength ap- 
proximation for the one-particle hamiltonian operator 

H = -ihv{(x x d x + a 2 d y ), (10) 

which acts on four-component wave functions of the form 

i> = $a-, ^b-) T , (11) 

where subscripts A and B correspond to two sublattices and subscripts + and — 
correspond to two inequivalent Fermi points. Thus, an effective long-wavelength 
description of charge carriers in graphene is written in terms of a continuum 
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model which is based on the Dirac-Weyl equation for massless electrons in 2 + 1- 
dimensional space-time, with the role of speed of light c played by Fermi velocity 
v « c/300 [HI EJ, see also Ref. [23]. 

In order to complete the Clifford algebra of anticommuting 4x4 matrices, 
one has to define 7 , 7 s , and 7 s = — Z7°7 1 7 2 7 3 (where 7 = 7°o:). It should be 
noted that hamiltonian (10) commutes with generators X& (k = 1,2,3) of the 
S , f/(2)-symmetry transformations, [T^, Tj\ = iekimT m , where 

T. = \l\ T 2 = i 7 5 , T 3 = i 7 3 7 5 , (12) 

and there is an arbitrariness in the choice of the representation of the Clifford 
algebra, which is due to a possibility of diagonalizing anyone of T^. A repre- 
sentation with diagonal 7 , in view of the block-diagonal form of a 1 and a 2 , 
see Eq.(9), corresponds to the choice of diagonal T 3 and may be denoted as the 
standard planar representation: 

7 " a 3 > 7 a 1 > 7 -a 2 



a 2 \ 5 . / -a 2 
a 2 J' 7 = \a 2 

Choosing T 2 to be diagonal, one gets the chiral planar representation: 
,3 J o 2 \ , ( -1 



(13) 



^ = -\l 2 J' 1 )■ ™ 

A rotation by angle i9 in the plane of a graphene sheet is implemented by 
operator exp(ii9£), where 

1 1 2 1 ( ^ 



E = — a V = - n 3 (15) 
2z 2 V -o- / 

is the pseudospin playing here the role of the operator of spin component which 
is orthogonal to the plane. The honeycomb lattice is invariant under the rotation 
by 2tt, 

exp (z27rS) if) = —if), (16) 
but is not invariant under the rotation by 7r, 

exp(i7rE)-?/> = 2iY<ijj, 
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i.e. under x — > —x and y — > —y. However, if the latter rotation is supple- 
mented by simultaneous exchange of sublattices, as well as Fermi points, then 
this combined transformation, 



is a symmetry one and can be regarded as the parity transformation for graphene. 
Note that transformation (17) is implemented by ia 3 = ry°7 3 in the standard 
planar representation, see Eq.(13), or by ry° in the chiral planar representation, 
see Eq.(14). The explicit form of R is extracted from Eq.(17): 



and it is given by —7 s in representation (13) or by — 7 3 7 5 in representation (14). 
Note that R is commuting with E (15) and H (10). 

3 Topological defects 

Topological defects in graphene are disclinations in the honeycomb lattice, result- 
ing from the substitution of a hexagon by, say, a pentagon or a heptagon; such a 
disclination warps the graphene sheet. More generally, a hexagon is substituted 
by a polygon with 6 — Nd sides, where Nd is an integer which is smaller than 
6. Polygons with N d > (N d < 0) induce locally positive (negative) curvature, 
whereas the graphene sheet is flat away from the defect, as is the conical surface 
away from the apex. In the case of nanocones with Nd > 0, the value of N d is 
related to apex angle 5, 

sm- = l- T , (19) 

and Nd counts the number of sectors of the value of ir/3 removed from the 
graphene sheet, see Fig. 3a. If Nd < 0, then —Nd counts the number of such sectors 
inserted into the graphene sheet. Certainly, polygonal defects with N d > 1 and 
N d < — 1 are mathematical abstractions, as are cones with a pointlike apex. In 
reality, the defects are smoothed, and N d > counts the number of the pentago- 
nal defects which are tightly clustered producing a conical shape; such nanocones 
were observed experimentally [24]. Theory predicts also an infinite series of the 
saddle-like nanocones with —N d counting the number of the heptagonal defects 
clustered in their central regions. However, as it was shown by using molecular- 
dynamics simulations [25], in the case of N d < —4, a surface with a polygonal 
defect is more stable than a similarly shaped surface containing a multiple num- 
ber of heptagons; a screw dislocation can be presented as the Nd — > — 00 limit of 
a 6 — iVd-gonal defect. 



i?exp(i7rS)^ = i {if) B -, 4>a-, 4>a+) T , 



(17) 




(18) 
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The two dimensional Dirac-Weyl hamiltonian on a curved surface with the 
squared length element ds 2 = ^y(r)dr J r- J takes form (see, e.g., Ref.[26j) 



H 



where 



and 



-ihvaHv) 



3 J '(r), (r) 



2g"\v)I 



(20) 
(21) 



r; fc (r) 



:Oc r. 



P ifc (r)a,(r 



2 glH ^ l d i9nk(r) + d k g nj {r) - d n g jk (r)} . 



In the case of a conical surface with a pointlike apex, one has 

1, g w = (1 - 77) V, 



(22) 



(23) 



where r and (p are polar coordinates centred at the apex, and — 00 < r\ < 1. The 
intrinsic curvature of a cone vanishes at r 7^ and possesses a 5 2 (r)-singularity at 
its apex (r = 0); parameter 77 enters the coefficient before this singularity term. 
Introducing ip' = (1 — rj)tp, one gets the metric in the (r, <//) coordinates, which 
is identical to that of a plane, but with ip' in the range < <p' < 2?r(l —77). Thus, 
quantity 2tti] for < t] < 1 is the deficit angle measuring the magnitude of the 
removed sector, and quantity —2tti] for —00 < i] < is the proficit angle mea- 
suring the magnitude of the inserted sector. In the case of graphitic nanocones, 
parameter r\ takes discrete values: 



T) = N d /6. 
Using Eqs.(21) and (23), one gets 

a r = a\ W = (1 - ^-V-V. 
It is straightforward to calculate the nonvanishing Christoffel symbols 



(24) 



(25) 



-pip 

rip 



-pip 

ipr 



and get the spin connection 

uj t = 0, uJip — —2(1 — 77)1]. 
Thus, hamiltonian (20) on a conical surface takes form 

H = -ihv {^dr + a 2 r- 1 [(1 - - zS] } . 



(26) 
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In the case of the planar graphene sheet (77 = 0), wave function (11) satisfyes 
antiperiodicity condition, see Eq.(16), 

iP(r, v? + 2tt) = -V>(r, </?), (27) 

i.e. the wave function is a section of a bundle with spin connection — 2E. 

Let us consider a graphene sheet with a pentagonal disclination (N d = 1). 
When circling once this defect, the two sublattices in the honeycomb structure 
are exchanged (see Fig. 3b), the two inequivalent Fermi points are exchanged as 
well. Circling twice this defect is analogous to circling once a hexagon in the 
honeycomb lattice without defects. The situation resembles that of a Mobius 
strip, where a double full turn is needed to arrive at the same point. Thus, in the 
continuum model description of graphene with a pentagonal disclination, wave 
function (11) has to obey the Mobius-strip-type condition: 



il>(r,<p + 2n)=iRil>(r,<p), (28) 

where R is given by Eq.(18), and, consequently, 

V>(r, ip + 4tt) = -ip(r, <p), (29) 

since R 2 = I. Note that the sign in the right hand side of Eq.(28) is chosen by 
convention. 

In a similar way, one can show that the wave function on a graphene sheet with 
a heptagonal disclination (Nd = —1) obeys the Mobius-strip-type condition as 
well. Moreover, it can be noted that sublattices (and Fermi points) are entwined 
in the case of odd N d and are left untwined in the case of even N d . Thus, the 
condition for the wave function on a graphene sheet with an arbitrary disclination 
takes form 

iP(r, <p + 2tt) = - exp (~i^N d R) ^(r, <p), (30) 

where the choice of sign in the exponential function agrees with the sign choice 
in Eq.(28). Our results remain unchanged if the opposite sign in Eq.(28) and, 
correspondingly, in the exponential function in Eq.(30) is chosen. 
By performing singular gauge transformation 

i>' = e in ij, n = V^R, (31) 
one gets the wave function obeying condition 

4>'(r,tp + 2n) = -,l/(r,(p), (32) 
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in the meantime, hamiltonian (26) is transformed to 



H' = e tn He- in = -ihv I a l d r + aV" 1 



(l-rjy^d^-i-rjR)-^ 



(33) 



where Eq.(24) is recalled. We conclude that a topological defect in graphene is 
presented by a pseudomagnetic vortex with flux Ndir/2 through the apex of a 
cone with deficit angle N d 7i/3. 



4 Ground state charge 

The density of states is defined as 

t(E) = -lm Tr(H -E-iO)~\ (34) 

7T 

where Tr is the trace of an integro-differential operator in the functional space: 
TrO = J <i 2 rtr(r|0|r); tr denotes the trace over spinor indices only. Since the 
continuum model of graphene without topological defects corresponds to the use 
of the free Dirac-Weyl hamiltonian in flat two dimensional space, see Eq.(lO), the 
density of states is immediately calculable and found to be proportional to the 
size of the graphene sheet 



S\E 
nh 2 v 

where S is its area. The ground state charge of the graphene sheet, 



r(E) - ^ (35) 



OO 

Q = ~\J dEr(E)sgn(E), (36) 

— oo 

is evidently zero, because Eq.(35) is even in energy. 

To consider the influence of topological defects in graphene on the density of 
states in the framework of the continuum model, we need the complete set of 
solutions to the Dirac-Weyl equation in this case 

H'ij>' = Eif>', (37) 

where H' is given by Eq.(33), and ip' is the spinor wave function obeying condition 
(32). In general, the contribution of a topological defect is added to Eq.(35), and, 
lacking the factor of area, it seems to be negligible. However, if this contribution 
contains a piece which is odd in energy, then the latter yields the nonzero ground 
state charge, see Eq.(36). In the following our interest will be in the search of 
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such a piece, and, as we shall see, its emergence is due to the appearance of an 
irregular solution to the Dirac-Weyl equation. 
Let us make unitary transformation 



i>" = U4>', H" = UH'U-\ 



where 



then 



u = u- 1 = —= 



URU- 1 = 



\_( I ia 2 
I 



-/ 

and transformed hamiltonian H" acquires a block-diagonal form: 



H" 



H x 
H-! 



(38) 
(39) 
(40) 

(41) 



where 



H s = hv <ia d r — a r 



' - //)" [isd v + ^rA+^a 3 



, s = ±l. (42) 



It should be emphasized that the definite sublattice (A or B) and Fermi-point (+ 
or — ) indices are assigned to the components of ip (11), while, after performing 
transformations (31) and (38), one gets ip" with components mixing up different 
sublattices and Fermi points. Certainly, the calculation of functional trace in 
Eq.(34) does not depend on the representation used, and it is just a matter of 
convenience to use a representation with the block-diagonal form of hamiltonian 
(41). 

Separating the radial and angular variables 



[r,ifi\E,n), 



(43) 



where 



(r,ip\E,n) 



(44) 



/ /n,i(r) e^ n+ ^ \ 

9n,l( r ) el{n+h)lp 

\ 9n,-i(r) e i{n -h) 

one rewrites the Dirac-Weyl equation as the system of equations for the radial 
functions 

Dl.\f f n , s (r) \ _ E ( f n ,s(r) 



J 







9n,s( r ) 



9n,s\r) 



(45) 
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where 



D n , s = hv[-d r + r 1 (l-r)) 1 (sn-rj)], 

Dl s = Hv[d r + r- 1 (l-rj)- 1 (sn+l-2rj)]. (46) 

A pair of linearly independent solutions to Eq.(45) is written in terms of the 
cylinder functions. In the case of | < rj < 1 (N^ = 3, 4, 5) the condition of 
regularity at the origin is equivalent to the condition of square integrability at 
this point, and this selects a physically reasonable solution. Thus, the situation 
is similar to that of rj = (absence of a defect), resulting in a density of states 
which is even in energy. In particular, it can be shown that the density of states 
in the case of rj — \ (Nd = 3) is given by Eq.(35), see Appendix A. 

In the case of < rj < \ (N d = 1, 2) and -\ < rj < (N d = -1, -2, -3) 
there is a mode, for which the condition of regularity at the origin is not equivalent 
to the condition of square integrability at this point: both linearly independent 
solutions for this mode are at once irregular and square integrable at the origin. 
To be more precise, let us define in this case 

n c = I [sgn(r?) - 1] , (47) 

and note that solutions to the Dirac-Weyl equation correspond to the continuous 
spectrum and, therefore, obey the orthonormality condition 

2w 00 

J dp J drr(l-r]){E,n\r,tp){r,ip\E',n') = 2 ^ ( ^__f\ w ,, (48) 


where a factor of 2 in the right hand side of the last relation is due to the existence 
of two inequivalent Fermi points. Then the complete set of solutions to Eq.(45) 
is chosen in the following form: 
regular modes with sn > sn c 



^ (r) ^ 1 ( Ji(i-r,)-i-F(kr) 1 l = s(n _ nc)) (4Q) 



V 9nA r ) J 2 v /tt(1 — 77) V sgn(£)J i(1 _„)-i +1 _ F (A; 
regular modes with sn < sn c 



fn,s( r ) 
9n,s( r ) 



- 1 ( fn^Mkr) \ _ 

2^(1-77) V sgn(E)J ni _ v) -i_ 1+F (kr) J' 

(50) 
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and an irregular mode 



9n c ,s(r) ) 2 v /tt(1 - rj) [1 + sm(2u E ) cos(Ftt)] 



/n e , s (r) \ = 1 



x 



X 



sin(z/£;) J_i?(/cr) + cos(^e) Ji?(/cr) 
sgn(E) [sin(v E ) Jx_ F (kr) - cos(v E )J_ 1+F (kr)} 



) 



(51) 



where k = \E\(Kv) 1 , ./^(u) is the Bessel function of order /z, and 



1 



-sgn(7/)+7? (I-77) 



(52) 
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Thus, the requirement of regularity for all modes is in contradiction with the 
requirement of completeness for these modes. The problem is to find a condition 
allowing for irregular at r — > behaviour of the mode with n = n c , i.e. to fix v E 
in Eq.(51). To solve this problem, first of all one has to recall the result of Ref. 
|27j . stating that for the partial Dirac hamiltonian to be essentially self-adjoint, it 
is necessary and sufficient that a non-square-integrable (at r — > 0) solution exist. 
Since such a solution does not exist in the case of n = n c , the appropriate partial 
hamiltonian is not essentially self-adjoint. The Weyl-von Neumann theory of self- 
adjoint operators (see, e.g., Ref. [2S|) is to be employed in order to consider a 
possibility of the self-adjoint extension for this operator. We show in Appendix 
B that the self-adjoint extension exists indeed, and the partial hamiltonian at 
n = n c is defined on the domain of functions obeying the condition 



where T(u) is the Euler gamma function, M is the parameter of the dimension of 
mass, and G is the self-adjoint extension parameter. Substituting the asymptotics 
of Eq.(51) at r — > (see, e.g., Ref. |29j) into Eq.(53), one gets the relation fixing 
parameter u E , 



Using the complete set of solutions, it is straightforward to determine the 
kernel of the resolvent, (r, <p\{H — oj)~ l \r' ', ip') (where u is a complex parameter 
with dimension of energy), calculate its functional trace, see, e.g., Ref. fl9j , and 
find density of states (34). Only the irregular mode contributes to the odd in 
energy piece of the density of states, which is given by expression 




(53) 




(54) 



t(E) 
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2(2F - 1) sin(FTr) 




2.F-1 / 

tan (f + | ) + (Jg, 


) COt (f + f ) 




7lE 


\( \E\ \ 
\Mv 2 J 


2(2F-1) / x 

tan 2 (f + s)-2cos(2F7r) + (Jg r ) 


2(1-2F) 

cot 2 (f + f ) 

(5 



Inserting Eq.(55) into Eq.(36), we calculate the ground state charge, 

Q = esgn [(l - 2F) cos 6], (56) 

where 

, , f sgn(w), n^O 
Sgn o^ = ( 0, u = 

The charge, if any, is accumulated in the vicinity of the defect, and its density is 
given by expression, see Ref. [PT] . 



p(r) 



2 sin(FTr) 



dww \K\-p(w) — Kp(w)] 



„»(! - ,),* 7 J^)— tan (| + J) + (^)^cot (f + J) 



(57) 



decreasing as an inverse power at large distances from the defect; here K^(u) is 
the Macdonald function of order //. 

In the case of rj < — ~ (iV d = —4, —5, . . .) there are two or more irregular 
modes, unless rj = — 1 (Na = —6). The case of more than one irregular modes 
will be considered elsewhere, while in the case of r\ = — 1 the irregular mode 
appears at n = n c with n c = —2s and has the form 



fn c ,s( r ) 



1 



sin(z/£) J_i (kr) + cos(^e) Ji (kr) 



9n c ,s( r ) J 2y/ir(l-ri) \ sgn(E) sm(u E ) Ji (kr) - cos(u E ) J_i (kr)} 

' (58 

hence this case corresponds to F = \ in Eq.(51), yielding the vanishing ground 
state charge. Our results are summarized in the Table. 

If the sign in the exponential function in condition (30) is changed to the oppo- 
site, then this corresponds to change F — > 1 — F. Our results remain unchanged, 
if, in addition, shift O — > O + n is performed. 



5 Discussion 

In the present paper we study the electronic properties of the carbon monolayer 
(graphene) with disclinations, i.e. 6 — A^-gonal (Nj ^ 0) defects inserted in the 
otherwise perfect two dimensional hexagonal lattice. The effects of the variation 
of the bond length or the mixing of tt- with u-orbitals caused by curvature of the 
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lattice surface are neglected, and our consideration, focusing on global aspects 
of coordination of carbon atoms, is based on the long-wavelength continuum 
model originating in the tight-binding approximation for the nearest neighbour 
interactions. Our general conclusion is that the dependence of the electronic 
properties on the value of Nd is not monotonic, but rather abruptly discontinuous. 
For some values of N d the density of states is predicted unambiguously by the 
theory, whereas, otherwise, its theoretical prediction involves some parameters 
which should be determined from the experiment. 

As it was already noted [20J , a defect with odd Nd entwines two sublattices, as 
well as two inequivalent Fermi points, and in the present paper we show that the 
correct condition for the electronic wave function involves operator R commuting 
with the hamiltonian, see Eqs.(30) and (18). However, much stronger impact 
on electronic properties might be drawn by the fact that for certain Nd irregu- 
lar modes emerge among the eigenmodes of the hamiltonian. It is instructive to 
compare two cases when the density of states remains almost the same as for pla- 
nar graphene, but for different reasons. In the case of the three-pentagon defect 
(Nd = 3), there is entwinement of sublattices and there is no irregular modes; the 
density of states is calculated (see Appendix A) and shown to coincide exactly 
with that of planar graphene. In the case of the two-pentagon defect (Nd = 2), 
there is no entwinement and there is an irregular mode; the density of states is 
even in energy (owing to F = ~, the odd piece vanishes, see Eq.(55)) and almost 
coincides with that of planar graphene. Thus, we disprove the controversial as- 
sertions that the density of states at the Fermi level is nonzero either at N d = 2 
[20] or at Nd = 3 [21]. 

The same unambiguous predictions are obtained for a graphene sheet with the 
two-heptagon defect (Nd = —2) and a graphene sheet with a dodecagon or six 
heptagons (Nd = —6): the density of states almost coincides with that of planar 
graphene. Evidently, the ground state charge is zero in all above cases. 

Let us turn now to the cases when our predictions are not unambiguous, since 
they involve self-adjoint extension parameter 0. These cases include graphene 
sheets with following defects: one pentagon (Nd = 1), one heptagon (N d = — 1), 
and three heptagons (Nd = —3). Contrary to the assertions in Refs. (20J, [21] , the 
density of states at the Fermi level is characterized by a divergent, rather than 
the cusp, behaviour in these cases, see Eq.(55): 



Actually, there are three possibilities: cos 6 > 0, cos 9 < 0, and cos 6 = 0. The 




-3. 



-1 



1 



(59) 
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question of which of the possibilities is realized has to be answered by experimen- 
tal measurements. First, the density of states in the vicinity of the Fermi level can 
be measured directly in scanning tunnel and transmission electron microscopy. 
Secondly, the ground state charge can be measured also, and our prediction, see 
Eq.(56) or the Table, is 

Q\N d =l = -Q\ Nd =-l = Q\N d =-3, (60) 

while the ground state charge density decreases by power law at large distances 
from the defect, see Eq.(57), 



p{r) 



9sin(7r/5) r (lj) r (lTi) I ( h \ f f ( 6 , tt \ at _-, 

20tt 2 r(|) r 2 KrMv) LVh \ 2 ^~ 4 ) ' ^d,— -*-, 



OO 



9sin(27r/7)r(t)r(il) ^(^) f tan(f + |), N d = -1, (61) 



— e — ^-y 



35tH rff 



V3 r (l) r (l) 1 ( h \\ + /0 I 7r\ O 



The results for the ground state charge in the case of cos > agree with 
the results of numerical atomistic calculation of the bond network with the use 
of recursion methods [30]. The pentagonal defect, as well as the three-heptagon 
one, is attractive, and the heptagonal defect is repulsive to electrons. The charge, 
negative or positive, is accumulated around the defect, and, at large distances 
from it, the decrease is the strongest one for a pentagon and the weakest one for 
three heptagons. 

It should be noted that at cos ^ and F ^ \ scale invariance is broken, and 
the appearance of parameter M with dimension of mass evinces this. In general, 
irregular mode (51) diverges at the origin as r~ v with v < 1. Scale invariance is 
respected by the condition of minimal irregularity 



e={ %™ 0i %\ !<£<f- (62) 
\ -f (mod2vr), | < F < 1, v ' 

which restricts the behaviour of the irregular mode at the origin to r~ v with 
v < |. Thus, both scale invariance and minimal irregularity favour definitely the 
choice of cos© = 0, when the density of states and the ground state charge are 
trivial. It would be inspiring, if the experiment could prefer other choices. 
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Appendix A 

In the case of r\ = | hamiltonian H s (42) takes form 



H, = hv 



( 3 1 
ia 2 d r - a 1 ^ 1 2isd v H h -<t : 



2 ' 2" J \ ' ^ A ' 1 ^ 
The kernel of the resolvent (the Green's function) of the hamiltonian is presented 



as 



(r, <p|(if a - 1 \r 



where the radial components satisfy equations 

-u hv(d r + r _1 2sra) \ f a l ^{r; r') af(r, r') 

hv [-d r + r'\2sn-l)} -u )\ a 12 (r, r') a 2 n 2 (r, r') 

—lo Hv(d r > +r'~ 1 2sn) \ / aj^(r, r') a* 2 (r, r') \ _ 

hv [-d r , +r'- 1 {2sn-l)] -u J \af (r, r') a 22 (r, r') 



(A3) 



note that a factor before the delta-function is due to (detg^y) 2 = [(1 — 77) r] 1 . 
The radial components can be found in the form 



dpp 



h 2 v 2 p 2 -u 2 ™' p 



(A.4) 



where 



a n,p( r , r ') 

22 / l\ 
a n,p{ r > r ) 



2u J2sn-l(pr)J2sn-l{pr'), 

2hvpJ 2s n(pr)J2sn-i(pr'), 

2hvp J2sn~l(pr)J2sn(pr'), 
2U) J2sn{pr)J2sn{pr'). 



(A.5) 
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Putting if' = ip and taking trace of matrix (A. 2), one gets 

tr(r,ip\{H s - w) -1 |r',¥>) = — [a"(r, r) + a^ 2 (r, r' 



oo 



= — / fc9 9 V" [^n-l(^) J 2 sn-l(pr') + J2sn(pr)J2sn(pr')] ■ (A. 6) 

The summation is performed with the use of the Neumann's addition theorem 
(see, e.g., Ref. [2H]), yielding the expression, 

oo 

tr(r^\(H s -u)-y^) = ^l h2 JP p _ - J \p(r - r% (A.7) 

o 

which is badly divergent at r' — > r. To tame the divergence, we define regularized 
kernel 



-X 



(r,<p\{H a -u) exp{-tH a )\r ,(p) = — ~ 



l£ \< p (r,r') a% p (r,r>))> ^.8) 

where t > is the regularization parameter. Now the limit r' — > r can be taken 
safely, yielding 



tr(r,<p\(H s -cu)- l exp(-tH 2 ) \r,<p) = ^J 





uj f dp p exp(—th 2 v 2 p 2s ' 



2nh 2 v 2 

where 



h 2 v 2 p 2 — U! 2 
£i(-^ 2 ), (A.9) 



w „ , . 2 



00 

*(«) = /-e- 



is the exponential integral (see Ref. [29J ) . Note that, actually, we have reiter- 
ated the derivation for the case of the planar graphene sheet (77 = 0): the only 
difference is that in the latter case all factors of 2 (including those at the order 
of Bessel functions) in the right hand sides of Eq.(A.5) are absent. 
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Since Eq.(A.9) is independent of r and ip, the integration over the surface 
yields a factor of its area: 



Tr(H s - c)- 1 exp (-tH 2 s ) = ^-^E 1 (-tu; 2 ). (A.10) 

The divergence of the last quantity in the limit t — ■> + does not contribute to the 
density of states, Eq.(34). This is due to a specific form of a discontinuity of the 
exponential integral at negative real values of its argument, lmE 1 (—u =F iO) = 
±in (u > 0). Consequently, we get finite result (35). 



Appendix B 

The partial hamiltonian corresponding to n = n c takes form, see Eqs.(45)-(47) 
and (52), 

h * ( d r + r~\\-F) \ ( . 

h = to{_ dr _ r -i F o J' (R1) 

Let h be defined on the domain of functions £°(r) that are regular at r = 0. Then 
its adjoint which is defined by relation 

oo oo 

I drr(l- V )[h^(r)n°(r) = J dr r(l - h?(r) (B.2) 

o o 

acts on the domain of functions £(r) that are not necessarily regular at r = 
0. So the question is whether the domain of definition of h can be extended, 
resulting in both the operator and its adjoint being defined on the same domain of 
functions. To answer this, one has to construct the eigenspaces of with complex 
eigenvalues. They are spanned by the linearly independent square-integrable 
solutions corresponding to the pair of purely imaginary eigenvalues, 

fct£±( r ) = ±iMv 2 £ ± {r), (B.3) 

where Mv 2 is inserted for the dimension reasons. It is straightforward to show 
that only one pair of such solutions exists 

£± (r) -±( e^K F (rMv/h) \ 

^ {r) ~ N V e^K^rMv/h) ) ' {BA) 

where N is a certain normalization factor. Thus, the deficiency index is equal to 
(1,1), and, according to the Weyl-von Neumann theory of self-adjoint operators 
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(see Ref. [28]), the self-adjoint extension of operator h is defined on the domain 
of functions of the form 

( Izt] ) = ^ 0(r) + c [e+(r) ~ e " er(r) ] ' (R5) 

where c is a complex constant and G is a real continuous parameter. Using the 
asymptotics of the Macdonald function at small values of its argument (see Ref. 
[29]). we get 



/ nc>s (r) \ = 2ce~ j T / -sin (f + f) 2 F T(F)(rMv/K)- F 
9n c ,s(r) J r ~* iN V cos (f + f ) 2 1 " i T(l - F)(rMv/H)- 1+F ) ' 

(B.6) 

which can be rewritten in the form of Eq.(53). 
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b— 



a) 




b) 

Figure 1: The planar honeycomb lattice as a composition of two triangular 
sublattices. The primitive cell with basis vectors Ci and c 2 is depicted in a), the 
triads Uj and Vj connecting different sublattices are depicted in b). 
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Figure 2: The first Brillouin zone is a hexagon with opposite sides identified, 
and, therefore, next to neighbouring corners are equivalent; two inequivalent ones 
can be chosen as lying on a vertical line. 
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b) 

Figure 3: Formation of a topological defect in graphene: a) one, two, or three 
sectors of 60° are removed from the lattice, b) if one sector is removed, then sites 
of different sublattices are identified, a and a', b and b', c and c', etc. 
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Table: 

The ground state charge in the case of the existence of the one irregular 
mode in the set of eigenmodes of the hamiltonian 



Nci 


2 


1 


-1 


-2 


-3 


-6 


F 


1/2 


1/5 


5/7 


1/2 


1/3 


1/2 


Q 
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